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It has been pointed out recently that the presence of dilaton field in the early Universe can dilute 
the neutralino dark matter (DM) abundance, if Universe is not radiation dominated at DM decou- 
pling, due to its dissipative-like coupling to DM. In this scenario two basic mechanisms compete, 
the modified Hubble expansion rate tending to increase the relic density and a dissipative force 
that tends to decrease it. The net effect can lead to an overall dramatic decrease of the predicted 
relic abundance, sometimes by amounts of the order of 0(1O 2 ) or so. This feature is rather generic, 
independent of any particular assumption on the underlying string dynamics, provided dilaton dom- 
inates at early eras after the end of inflation but before Big Bang Nucleosynthesis (BBN) . The latter 
ensures that BBN is not upset by the presence of the dilaton. In this paper, within the context 
of such a scenario, we study the phenomenology of the constrained minimal supersymmetric model 
(CMSSM) by taking into account all recent experimental constraints, including those from the LHC 
searches. We find that the allowed parameter space is greatly enlarged and includes regions that 
are beyond the reach of LHC. The allowed regions are compatible with Direct Dark Matter searches 
since the small neutralino annihilation rates, that are now in accord with the cosmological data on 
the relic density, imply small neutralino-nucleon cross sections below the sensitivities of the Direct 
Dark Matter experiments. It is also important that the new cosmologically accepted regions are 
compatible with Higgs boson masses larger than 120 GeV, as it is indicated from the LHC exper- 
imental data. The smaller annihilation cross sections needed to explain WMAP data require that 
the detector performances of current and planned indirect DM search experiments through 7— rays 
should be greatly improved in order to probe the CMSSM regions. 
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I. INTRODUCTION 



Cosmological data accumulated from various observations over the past twelve years leave little doubt 
that Dark Matter (DM) and Dark energy (DE) occupy the major portion of the total matter-energy 
budget of the Universe [iKa] . These results follow from best-fit analyses of various astrophysical data to 
the Standard Cosmological Model (ACDM) which can successfully describe the evolution of our Universe. 
The model is based on a Friedmann-Robertson- Walker (FRW) cosmology, which involves cold DM, at a 
percentage 23%, baryonic matter at 4% and a positive cosmological constant A > which is introduced 
in an ad-hoc manner in an attempt to describe the vacuum energy density. 

The WMAP data @, 0| , with their unprecedented accuracy, have constrained the predictions of various 
particle physics models, among these supesymmetric (SUSY) theories that predict the existence of WIMP 
particles, the neutralinos, that would play the role of DM. Besides on-going LHC experiments running 
at y/s — 7 TeV, with luminosities that are constantly improving, have constrained the Standard Model 
(SM) Higgs mass to a narrow window while they have put new limits on the sparticle masses of SUSY 
theories [6j. The latest CERN LHC data could be interpreted as a first evidence for the existence of 
a Higgs with a mass around 125 GeV 0- This may have severe implications for all SUSY models, 
especially if combined with the WMAP data. 

After one year operation the data provided by the LHC give us the opportunity to delineate regions 
of the SUSY parameter space that are also compatible with the cosmological data. In particular, the 
leading supersymmetric DM candidates, the neutralino 8], still lacking experimental verification, results 
to a DM abundance which is severely restricted by cosmic microwave background data and various data 
from collider experiments, if calculated in the context of the simplest supersymmetry models (minimal 
supersymmetric model embedded in minimal supergravity [9] ) . In calculating neutralino DM abundances 
it is usually assumed that radiation dominates during DM decoupling. However, there are various mech- 
anisms originating from string theories, due to dilaton dynamics, which assume that our Universe is not 
radiation dominated during DM decoupling. These may alter the predictions for the DM abundance in 
a dramatic way. 

In general the existence of scalar fields in the primordial Universe, that contribute to the energy density, 
may play a significant role and upset the conventional picture. Among these are quintessence scenarios 
[lOj that have been invoked in an attempt to explain the vacuum energy, whose existence affects relic 
abundances since Universe is not radiation dominated during DM decoupling. In these models DM relic 
density is enhanced [ll|, and in some cases this enhancement can reach the level of ~ 10 6 or so fl2T ]. 
For a review see for instance (l3| . In general, modifications of the expansion rate which deviate from 
the standard cosmological scenarios may have dramatic consequences for the DM relic density [3, [l5[ • 
Inversely, it has been pointed out that the observed amount of DM puts rather strong constraints on 
possible modifications of the Universe expansion rate at early eras [la] , A particular subclass of these 
models, the trackin g qu intessence scenario, assumes that the quintessence field is in a kination-dominated 
phase at early eras [17| • In this context the predictions for the gravitino and axino DM are considered in 
[l8| while in (l9l | the predictions for the neutralino DM relic, in the popular supersymmetric schemes, is 



discussed in the light of the observed e -spectrum by PAMELA [20( and Fermi-LAT observations [21 1 



Also in other string inspired models with a time-dependent dilaton sources 22], whose evolution obeys 
a non-equilibrium string dynamics 23[, the amount of thermal neutralino relic abundance is diluted by 



factor of O(10), relative to that calculated within the A CDM-minimal supergravity cosmology and such 
models seem to survive the stringent tests of LHC [24| . In this case the dilution is due to the appearance 
of a friction-like term, on the right-hand side of the appropriate Boltzmann equation. This term plays 
also a significant role in other considerations studied in [25j . 

In this paper we explore the phenomenological consequences of the constrained minimal supersymmetric 
standard model (CMSSM) model, which we consider as a prototype, in order to study the effects of the 
dilaton field assuming its energy prevails over radiation at early eras. In ref. |2q ] it was shown that the 
mechanism for the dilution of DM relic abundances is more general relying on more generic features of the 
dilaton dynamics which are independent of the non-criticality [13] of the underlying string theory. For 
dcfiniteness we model the dilaton evolution at early eras to be dictated by exponential-type potentials, 
which arc rather generic in string theories and arise from loop corrections, occurring also in a wide class 
of quintessence scenarios, inflation models and supergravity theories. The presence of such a dilaton field, 
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that dominates over radiation long before nucleosynthesis, affects relic abundances in a dramatic way due 
to its dissipative like coupling to Dark Matter. Just before and during BBN era dilaton's energy is much 
smaller than radiation energy while dilaton attains a constant value at hadronization prohibiting such 
a dissipative coupling to ordinary hadrons. Therefore, dilaton has no effect when hadrons are formed 
and abundances of ordinary matter are not affected, unlike DM relic abundances which can be smaller 
by factors as small as 0(1O~ 2 ) , as compared to the conventional calculations. This allows for smaller 
annihilation cross sections, in the popular supersymmetric schemes employed in literature, and enlarges 
the cosmologically allowed parameter space, which is moved to regions that would be otherwise forbidden, 
altering the potential of discovering supersymmetry at collider experiments. We should remark that in 
this scenario the cross sections of any processes are not affected by the presence of the dilaton, it is the 
relic density that it does due to the modification of the Boltzmann equation. As we shall see this causes a 
dilution, in general, and therefore agreement with WMAP data is obtained for smaller annihilation cross 
sections. This in turn implies smaller elastic neutralino-nucleon cross sections which are of relevance 
to other experiments aiming at discovering DM (for reviews see [HI), like direct [2!| [3(| and indirect 

mm 13 dm searches ( f ° r 

a review see [32|). It is known that direct DM experiments impose constraints 
on supersymmetric models by excluding part of their parametric space that is close to the focus point 
region. Therefore a phenomenological study requires that both WMAP bounds and those stemming from 
direct DM searches are considered. 

In this paper modeling the dilaton evolution at early eras as prescribed in [26j | , we follow a phenomeno- 
logical analysis by taking into account the recent accelerator data from LHC, the constraints provided 
by WMAP7, and DM direct searches. The current data from indirect detection of DM through 7-rays 
are briefly discussed. In this discussion we do not consider energy spectra of fluxes from positrons and 
anti-protons. Although these can be conditionally be interpreted as due to annihilation of DM [33T - l35j 
nevertheless more conventional explanations exist [36]. In fact, the former may be explained as emission 
from pulsars and the latter are compatible with standard production mechanisms from Cosmic Ray which 
impinges the stellar gas. 

The paper is organized as follows: 
In section II we outline the salient features pertinent to the dilaton evolution and its dominance at early 
eras while in section III we discuss the role of the dilaton field and its effect on the DM relic abundance 
paving the ground for the phenomenological study of the constrained MSSM model (CMSSM) which is 
done in section IV. In section V we discuss the importance of 7-ray indirect DM searches for this scenario. 
Our conclusions are presented in section VI. 



II. EARLY DILATON EVOLUTION AND ITS DOMINANCE 



Since the main hypothesis is dominance of the dilaton energy, over that of radiation and matter, one 
can omit their corresponding contributions to the equations of motion. Therefore, the equations that 
govern the evolution of the dilaton field take on the form 

4> + 3H4> + v'(cj)) = 0, 

3H 2 = ^ + V&), 

2H = -(£> +P ) - -^ 2 . (1) 

In these the field (f> is dimensionless and the potential carries dimension mass 2 . The first of these equations 
is not independent but is derived from the other two. We further assume that the evolution of the dilaton is 
dictated by an exponential-type potential having the form V ~ k ^ . Such potentials are motivated by 
inflation (37j , quintessence scenarios (38j and they can also occur in string theories as perturbative or non- 
perturbative corrections. Such potentials accept dilaton solutions which are logarithmically dependent 
on the cosmic scale factor a(t) , during early eras, having the form 

4> = c In (£\ + (f)! . (2) 
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In this c is a constant and aj = a(tj) is the cosmic scale factor at the maximal reheating temperature, 
which was reached after inflation, denoted hereafter by Tj. The corresponding time at which this occurred 
is tj. 

Evidently the above solution for the dilaton holds in epochs tj < t < tx or, in terms of the cosmic 
scale factor, when aj > a > ax- Throughout the subscript X signals the end of the period in which 
dilaton dominates over radiation and matter. At tx dilaton stops evolving, it attains a constant value 
so that its kinetic energy vanishes, and its potential energy becomes very small so that radiation energy 
starts taking over. 

Both ai and ax are inputs and it is more convenient to trade them for the parameters r and b dchncd by 
r = ln(a//oo) and b = ln(ax/ao) respectively. The first is set by the choice of the reheating temperature 
Tj . In a supersymmetric model with the MSSM content, whose sparticle mass spectrum is in the range of 
a few TeV, supersymmetric as well as SM particles are all relativistic and the effective number of degrees 
of freedom is g e ff = 228.75, independently of the precise sparticle mass spectrum. This entails to a 
value r — —50.86 — ln(Tj/10 9 GeV). Throughout we assume that Tj is of the order of 10 9 GeV and 
hence r is in the vicinity of ~ —50. Although our conclusions are not sensitive to the value of Tj, as 
long as it is much higher than the TeV scale, i.e. the typical supersymmetry breaking scale, in this work 
for definiteness we consider values of it in the 10 9 GeV range which are supported by CMB constraints, 
which put lower bounds on the inflationary reheating temperature (3£^, of the order of ~ 10 6 GeV, and 
Baryogenesis through thermal Leptogenesis which demands temperatures ~ 10 9 GeV |40j . 

As for the other parameter b, a first upper bound arises from BBN which should not be affected 
by the presence of the dilaton. Therefore, the end of the dilaton dominance period tx must occur 
before Nucleosynthesis. In particular just after tx radiation starts contributing and overwhelms dilaton's 
energy at the time of Nucleosynthesis. The latter took place at Tbbn — 1 MeV, corresponding to 
\il{<ibbn /oo) — —22.5 and therefore b must be less than this. However the bound put on b must be 
even smaller since dilaton should get an almost constant value when hadrons are non-relativistic and this 
occured before Nucleosynthesis. Its constancy is rather mandatory in this era, otherwise the diluting 
dilaton mechanism will affect the abundances of the known hadrons and especially nucleons which we 
do not want to occur. In fact, the couplings of dilaton to matter density appear through dissipative 
terms ~ (g m — 3p m )0, which modify the corresponding continuity equation for matter. Such terms 
are vanishing when hadrons are relativistic, that is at temperatures higher than about Th ~ 1 GeV, 
corresponding to ln^^/ao) ~ —30, causing no harm. However below Th hadrons are pressureless and 
dilaton couples to hadrons through ~ g m cb. Therefore, in this temperature regime the dilaton has to be 
almost constant in order to suppress its coupling to hadronic matter. A reasonable value is Th = h-QCD, 
with Aqcd — 260 MeV the characteristic QCD scale. This lowers ax, and hence b, suggesting a value 
given by b ~ —28.4. The interesting point is that for such values of Th the coupling of the dilaton to 
supersymmetric matter is non- vanishing during DM decoupling. In particular, if DM has supersymmetric 
origin, its decoupling occurs at a temperature between Tdm — 5 — 20 GeV, corresponding to values of 
the cosmic scale factor ln(aj)M /clq) in the range ~ —31.5 to ~ —33.0, when dilaton is still the dominant 
source of energy. This may have dramatic effects for the DM relic abundances as we shall see. 

As a side remark, note that a constant dilaton in the range t > tx cannot account for changes of the 
fine structure constant, Aa/a ~ 10 -5 , over cosmological time scales. The constancy of the dilaton in this 
era naturally implies small quantum fluctuations A (f> <C 1 which may put limits on the dilaton-matter 
coupling approaching the capability of Eotvos-like experiments. 

A dilaton solution given by Eq. ((5J) implies that the time derivative of <j> is linearly related to the 
Hubble rate through cj) = cH , for t < tx- Then by the third of equations {I), which can be solved for 
H , we find that H is inverse proportional of the time t. Solving H = a/a we can get the expansion rate 
a{t) and from this the form of the dilaton field as function of the time, 

In this subscripts I denote quantities evaluated at tj and t„ =t — tj. 

Knowing the dilaton and the Hubble rate, from the second of Eqs. (JTJ) the dependence of the potential 
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on the values of the expansion rate and the dilaton at the reheating time can be derived 



V(<f>) = I Hfe-'W. (4) 

Using this one finds that the value of the potential at ax is exponentially suppressed relative to that at 
a/. In fact 

V{4>x)/V{<t>i) = exp[-c(<f> x - <pi)} -exp^Bc 2 ). (5) 

In this the parameter B is B = b — r, with &, r defined before, and therefore B has a value around 
~ 20. For a positive potential, which drops as Universe expands, the constant c should be bounded by 
c 2 < 6 , as is evident from ((4]), while as we shall see, dominance of the dilaton energy over radiation is 
achieved for c 2 > 4 . Since B has a value ~ 20, by Eq. J5]) we see that dilaton's potential energy is 
exponentially suppressed long before Nucleosynthesis. In fact, the ratio of the potential energy, at the 
end of the dilaton-dominated era, to the same energy at reheating temperature, given above, is found to 
drop by at least fourty orders of magnitude in this scenario, and the same holds for the total dilaton's 
energy since the ratio of its kinetic to its potential energy is constant. This is easily derived by using 
the second of Eq. ([IJ , using the fact that <j) = cH . Hence in this model the barotropic index w has 
a constant value, in the regime tj < t < tx , given by w — c 2 /3 — 1 . From the limits put on c 2 we 
conclude that w lies in the range 1 < w < 2 . 

For the dilaton to radiation energy density ratio one finds that 



(>r \Pr 



(G) 



where {p(f,/pr)\i is the same ratio at the reheating temperature Tj . In Eq. ([5]) we have reinstated 
dimensions and hatted densities carry dimension energy 4 . This ratio ought to be much larger than unity 
in the entire range tj < t < tx , since dilaton dominates over radiation in this period, and should decrease 
as Universe expands so that eventually before Nucleosynthesis radiation starts dominating. For this to 
happen the constant c should satisfy c 2 > 4 . Combined with the bound c 2 < 6 discussed earlier, which 
ensures that we have a positive potential that drops as Universe cools, we see that c should lie in the 
rather tight range 4 < c 2 < 6 . As a remark we point out that the ratio {p<j>/p r )\i turns out to be 
proportional to ( Hi / Hq) 2 , where Hj is the value of the Hubble rate at the end of inflation and Hq 
its corresponding value today. Measurements of the power spectrum of scalar and tensor perturbations 
yield bounds on the inflationary potential which in turn imply an upper bound for (Hj/Hq), (4lj . 
However this bound is rather weak in imposing a severe upper bound on the ratio {pip/ p r )\i which can 
be comfortably much larger than unity, as demanded, conforming at the same time with the constraints 
arising from the power spectrum measurements (for details see [26(). 



III. DILUTION OF RELIC DENSITIES 



Due to the presence of the dissipative term the continuity equation for the energy-matter density 
includes an additional term, 

^ + 3H(p + p)- A (p _ 3p )=0. (7) 

In this equation the division by y/2 of the last term is due to the normalization of the dilaton whose 
kinetic energy in our notation appears as p^ zn = (fi 2 /2 as is seen from the second of Eq. ([TJ. It should 
be noted that throughout this work we are assuming the lowest order, in a' , contributions to the form 
factors e~^^ and Z(4>) associated with the scalar curvature R and dilaton kinetic terms of the string 
effective action, in the string frame, and hence the simple expressions for the ^-dependent term of Eq. 
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([7]). Also for simplicity the dilatonic charge has been taken vanishing. We are aware of the fact that 
the couplings of the dilaton to matter may evolve in time with the dilaton itself and depend on the 
particle species in a non- universal way. Therefore, other options are available [42], [43[. Besides there may 
also exist additional contributions to the continuity equation which depend on the dilatonic charges of 
the species under consideration as studied in |42h44|. Certainly in order to take into account all these 
effects one needs a better understanding of the underlying string dynamics including perturbative and 
non-perturbative string effects which are lacking. Therefore, in this work we consider the simple case 
scenario in which the gravi-dilaton effective action is approximated by its lowest order Lagrangian in the 
string slope a' . 

The last term in Eq. (J7J drops when matter is relativistic. Obviously this is also the case for radiation 
since the difference p — 3p stays zero during the entire cosmological evolution. Therefore, for particles of 
a certain mass m the last term drops for temperatures T 3> m and the effect of the dilaton dissipative 
force is absent. However, as temperature drops this term eventually starts contributing for all period 
until the dilaton attains its constancy which occurs after DM decoupling but before hadron formation, 
and thus before BBN, in the way described in the previous section. Including the collision term and 
neglecting the pressure, for temperatures in the regime T < m , Eq. ([7J implies the following equation 
for the number density 

^ + 3Hn + (av)(n 2 -n 2 eq )- ^=n = 0. (8) 

This is suitable for describing the evolution of the number density during periods for which the particle 
is non-relativistic and pressure practically vanishes. During eras in which the particle is relativistic the 
last term in Eq. ((5J) is absent, as we argued before, since the last term in the continuity equation (J7J) is 
absent too. Thus Eq. ((8} receives the well-known form of Boltzmann equation [45| . Using the number to 
entropy density ratio, Y = n/s , this equation takes on the following form 

1/2 

m(av) ^—JLg eff j (h+~—)(Y 2 -Y 2 q ) + S(x) Y, (9) 

where x stands for x — T/m , with T the photon gas temperature which is related to the radiation 
density p r through 

Pr = 30 5e//( ^ ^ ' 

The function g e ff(T), at a given epoch, counts the degrees of freedom of the particles that are relativistic 
that epoch. In Eq. ([9]) Gjv is Newton's constant and the quantity h stands for the cntropic degrees of 
freedom related to the entropy density through s — 2tt 2 T 3 h(T) / 45 . Note the appearance of the the 
prefactor £(x) in the first term of Eq ^ which is given by 

\ p r SnGNPrJ 

In this p r ,Pm and p^ are the radiation, matter and dilaton energy-densities respectively. Recall that we 
use a dilaton density having dimensions m 2 , see Eq. [TJ Note also that no cosmological-term contributes 
to since such a term is absent during the early Universe evolution which includes the DM decoupling 
era. £(x) is smaller than unity having therefore the tendency to decrease relic abundances at slower rates 
as compared to the conventional cases, where this factor is unity due to the fact that in these approaches 
Universe is assumed radiation dominated during DM decoupling. Therefore, if it were only for £(x) 
an increase of the DM relic density, relative to the conventional scenarios, would have been predicted. 
However in the same equation an additional source term S(x) is present, given by 
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q(x) = n/T s /h 
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FIG. 1: The LSP Dark Matter number density to entropy density ratio q — as function of x = — ^ — in a 
particular supergavity model. The values of £, S denote the status of the ^-factor and the source respectively ( / 
for open, x for switched-off ). For comparison the corresponding equilibrium density qo{x) has been also drawn. 
(From Ref. 

In this expression <j)' denotes differentiation with respect ln(a/ao). If dilaton decreases, with decreasing 
temperature, then the source term is positive and acts in the opposite direction tending to decrease Y 
after DM decoupling. Therefore, two different mechanisms compete in this scenario and the net effect 
may be a decrease of the relic abundance. 

As an example, in Fig. Q] we display the effects of the the presence of the factors £(x) and S(x), as 
given by Eqs (fTOf and (fTTj). for the ratio of the number to entropy density n/(T 3 h). The displayed 
figure correspond to a supergravity model with inputs given by mo = 1100 GeV, M1/2 = 1200 GeV and 
Aq = GeV with tan/3 = 40 and the parameter /i > 0. For the particular SUSY inputs the LSP Bino 
has a mass ttilsp = 527.2 GeV. Eq §§§ has been integrated numerically which yields more reliable results 
than other approximate schemes employed in other works. The value of b = ln(ax/tto) , setting the onset 
of the epoch after which the dilaton is constant, has been taken —28.4 corresponding to a temperature 
Aqcd = 260 GeV as we have already discussed. For comparison in Fig. Q] except the ordinary case 
scenario, where both the source S and the £— factor are absent (red solid line), the cases where both 
£ and S are open (green dashed-dotted line), or when only the £— factor is present (blue short-dashed 
line) are also shown. The very thin dashed line, that rapidly drops, is the equilibrium density. In the 
case, that both terms are switched on, the density is monotonically decreasing after decoupling due to 
the appearance of the source term. The rapid change around x ~ 0.0005, corresponds to a value of the 
cosmic scale factor specified by b = —28.4 where dilaton reaches its constancy. In the specific example 
shown the relic density is diluted by a factor of ~ 50, as can be seen by comparing today's density 
predictions for the conventional case (red solid line) and the case where both £— factor and the source 
term are switched on (green dashed line). In the first case the relic density predicted is Olsp = 6.059 
while in the second case the relic density is considerably reduced falling into the WMAP allowed range 
^lsp/io = 0.1116. 

Although in this work we solve © numerically, which yields the most accurate and reliable results, we 
shall proceed to a semi-analytic treatment of © and illuminate some of its qualitative features that are 
worth pointing out. The presence of in this equation it effectively decreases the thermal average 
(a v) since £(x) < 1. As a consequence this term has the tendency to shift the freeze-out point Xf = Tf/m 
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to higher values, as compared to the conventional case where this term is unity, while the source term 
has little effect on Xf. We shall quantify this later. This results to a larger value of 7 at Xf. Shortly 
after freeze-out (av) n -C H and the first term in (JUJ) can be dropped. Then we can integrate from Xf 
to today's temperature and find 

= r s{x)dx ■ (i2) 

Integrations of the source S(x) can be performed analytically, if the entropy contributions in fjl If) are 
neglected, which is a decent approximation. In particular 

X2 S(x)dx= ^-f^ = iW|M2 =^(ln^ + iln^4p)), (13) 
V2 V2 V2 V T 2 4 g efs {T 2 )J 

where the particular dilaton solution ([2]) has been employed. If the lower limit x\ of the integration is 
below xx = Tx/m then T\ should be replaced by Tx on the r.h.s. Therefore from (fT2")l and (|T^|) we get 

InY(x) = lnY(x f ) - -£= fin f + \ In MM\ , (14) 

V2 V 4 9eff{xf)J 

for any x < Xf. Since c < the logarithmic term ~ In a; explains the linear drop of \nY(x) below Xf in 
Fig. Q] (green line). From this it can be seen that the value of Y reached at today's temperature is rather 
sensitive to the choice of the parameter b. In particular for two different values b and b' the resulting 
values are related by 

Y\x Q ) = Y(x ) exp C(b ^_ b) . (15) 

As the value of c is confined in the range 2 < \c\ < 2.45 we observe that by increasing b by a unit, Yq 
decreases relative to Y by factor of about of ~ 5. 

One can attempt to find semi-analytic solutions of Eq. @ which constitute approximations of the 
numerical results. Using standard techniques, already encountered in the dilaton-free case, where the 
terms £, S are absent one can solve (|9|) in a semi analytic way. This treatment leads to the following 
result for the relic density 1 

nhg = R(nh 2 ) no _ maton , (i6) 

where the factor R on the r.h.s. accounts for a change of the result of the conventional case. In this 
approximation R is given by 



R = k£ 1 {xf) exp 



S(x) dx 

XQ 



(17) 



where A: is a constant, with k < 1, which will be discussed below. Wherever appears Xf stands for the 
freeze-out temperature and xq is the value of x = T/m today. The integrations of the source S(x) in 
this equation can be performed analytically, as we have already discussed, using (fH?)) . 
In the approximate scheme that leads to (|16j) the freeze-out point is found by solving 

xj 1 = In (o.076 Cf(cf + 2)m M Planck x 1 / 2 (a v) f 9~ff 2 (xfj) 

+ In &f) - In (1 - 3x f /2 ~ xj (1 + c f ) S(x f )) , (18) 



1 Such a treatment has been already considered in the first of \'2'M . Here we refine and modify the results presented in that 
work and adapt them to the model at hand. 
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where Mp^^ = 1.22 x 10 19 GeV. When £ = 1 and 5 = the constant k becomes unity and we recover 
the standard result. In this equation the constant ct is of order unity and it is defined by Y — Y eq = cfY eq . 
In the standard cases the value of Cf is adjusted so that agreement with the approximate results is 
obtained. The source term at Xf is S(xj) = —c/y2xf and thus the last term in (fl~8f can be shown 
to be a small number if Xf is O(0. 01). The quantity £(x) is the inverse square root of the ratio of the 
dilaton to radiation energy at x, when dilaton energy dominates. Therefore the In £(2;/) term in (|18p 
is negative and it dominates over the last term in the same equation. Hence, the freeze out point is 
shifted towards larger values as compared to the standard result. For instance, with values of the dilaton 
to radiation energy ratio ~ 10 2 — 10 4 , in the region T = 5 — 50 GeV, one can find from (jT5J) that the 
relevant shift is at most 20%. This we have also verified numerically. Therefore for a WIMP, values in 
the range ~ 0.05 — 06 are expected for the freeze-out point, provided the ratio of the dilaton to radiation 
energy is not exceedingly large. 

Regarding the factor R of Eq. (TTT)) . the factor £ _1 (x/), and the constant k, arose by approximating 
the integral J Q f (crv) dx, encountered in integrating © from xq to xj, by £(x/) f^* (crv) dx. The 
constant k is actually the ratio of the approximated to the exact integral. Since £(x) > £(£/), for any 
x < Xf, by its definition the constant k is less than unity. Our numerical procedure shows that k takes 
values in the whole region from 0.1 to 1.0, with no preference to particular values or region, and it 
depends on the inputs in a way that it does not allow us to express it by a simple empirical formula. 
Hence R as it stands cannot be trusted to derive precise results and this is the main reason we decide 
to solve Boltzmann equation numerically. However R bears the main qualitative features, especially the 
exponential factor suppressing the relic density and it serves as an order of magnitude estimate of the 
actual result. In R, the factor £,~ 1 (xf) tends to increase relic density, as is evident from Eq. (TIT))) . The 
enhancement caused by £ _1 (x/) is followed by the exponential term in Eq. (1171) that acts in the opposite 
direction tending to decrease the relic density. 

R depends on Tj , Tx , which sets the end of the dilaton dominance era, and the value of the slope c, 
which are free parameters. Besides the ratio of the dilaton to radiation energy density A = p^/RttGn Pr 
at a given temperature, say Tx, which is shall denote by Axj is also a free parameter. Given Ax, by 
Eq. ([5]), its value at any other temperature can be derived. Although Tj is an arbitrary parameter the 
preferred values of it are around 10 9 GeV as dictated by inflation and Leptogenesis scenarios as we have 
already pointed out. The constant c lies in the narrow range 2 < |c| < 2.45, as we have already stated 
and Tx should take values larger than the typical hadronic scale, Th — 1 GeV or less, but not much 
less. In our numerical studies we use values Tx — ^qcd which yield the maximal suppression to the 
relic density. 

It facilitates discussion if R is cast in the following form 

R = ker p . (19) 
In this the exponent P, to a good approximation, is given by 

P = \{±-c 2 ~V2c) hag - iln(l + A x ). (20) 

From this we see that suppression of the relic density is obtained if P > 0. In P the prefactor of In ^ 
is positive in the allowed range 2 < |c| < 2.45, taking values varying from 1.41 to 0.73. Since Tj > Tx 
the first term suppresses R, and hence relic density, while the second term tends to increase it. The 
maximal suppression effect is implemented for the smallest possible values of |c|, Tx and Ax- As a case 
leading to a severe reduction of the relic density, if we take for instance Tx — 0.5 GeV, Tj ~ 20 GeV 
\c\ — 2.05 and Ax — 100 from the exponential term alone we get R ~ 0.065 which is of order 0(1O -2 ). 
The prefactor k decreases further this value. 

In the discussion that follows we shall present results for the relic density for a WIMP particle of 
mass m x derived by solving Boltzmann equation numerically, and compare it to that obtained by using 
the semi-analytic method based on Eqs. (|16I17[) . Although complete agreement between numerical and 
approximate results cannot be obtained, the comparison will serve in order to check the real magnitude 
of the reduction and compare it to that implemented by the factor R, which was discussed before. 
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a b 
[cm 3 /s] 


m x 
[ GeV] 


in 2 w r> n;, 2 l n;. 2 I 'QHq \num 

X-f 1U X It \num \LiIq \app D£,2 1 

^ ^*^0 a PP 


2.5 x 10~ 25 


500 
900 
5000 


0.0329 2.987 4.741 xlO" 2 4.062 xlO" 2 1.167 
0.0322 1.343 2.154 xl0~ 2 1.822 xlO~ 2 1.182 
0.0304 0.129 0.0020 xlO -2 0.0017 xlO -2 1.176 


1.0 x 10" 27 


250 
800 
5000 


0.0441 5.240 0.1542 0.1346 1.146 
0.0418 1.096 3.380 xlO -2 2.859 xlO" 2 1.182 
0.0387 0.092 0.273 xl0~ 2 0.230 xlO" 2 1.187 


1.0 x 10~ 29 


250 
800 
5000 


0.0592 3.465 7.656 6.678 1.146 
0.0552 0.741 1.711 1.453 1.178 
0.0500 0.064 0.143 0.122 1.172 


1.0 x 10 -29 1.0 x 10~ 29 


250 
800 
5000 


0.0589 3.490 7.345 6.574 1.117 
0.0549 0.746 1.674 1.432 1.170 
0.0498 0.064 0.138 0.119 1.159 



TABLE I: In the first two columns are the values of the coefficients a, b, in cm 3 /s, for the thermal 
averaged cross-section (a v) = a + bx while in the third column is the WIMP mass. In the next four 
columns we display the freezing point ( shifted value ) x c j, the reduction factor 10 2 x i?, the relic density 
derived by solving Boltzmann equation numerically, \num, and its approximate value fift 2 , \ app . In the 
last column their corresponding ratio ilh^ \ num / Qhy \ app is shown. 



Furthermore such a comparison may suggest possible refinements which will improve the approximation 
bridging the discrepancy between approximate and numerical results. 

The factor R is dominated by the exponent P, given by Eq. (|20j) . which however is approximate. Its 
exact form, as read from Eq. (fl~T|). is 



1 f Xf 
Pexact = --ln(l + A/)+ / S(x)dx 



= -I In 
2 



l + A x exp( C 2 -4) (fc + ln^ + iln 

V T 4 9eff{T )J 



In passing to the second equation the ratio of the dilaton to radiation energy at Tf, denoted by A/, 
has been expressed in terms of the corresponding ratio at Tx, named Ax, using Eq. ((6]). Also Tx has 
been expressed in terms of the parameter b = ln(a//ao). In this way Pexact is expressed in terms of the 
parameters b , c and Ax- The freeze-out temperature Tf appearing in Eq. (|21l) is read from (|18p . The 
temperature Tx appearing in the last line of this equation cannot be expressed in an analytic manner in 
terms of 6, since b = \n(T /T x ) + 1/4 \n(g e ff(T x )/g e ff{To))- However for T x in the range from A-qcd 
to 1 GeV, which we are mainly interested in, the quantity 1/4 ]n(g e ff(Tx)/geff(To)) varies little, taking 
values between 0.66 and 0.77, and one can replace it by its average ~ 0.71 reproducing satisfactorily the 
exact result. In order to fully quantify the factor R, and present an expression for it in terms of the 
parameters describing the model, the prefactor k appearing in Eq. (|17j) is also needed. This has been 
discussed before, in this section, and cannot be expressed in an analytic way. Later we shall give an 
approximate formula which reproduces fairly well the results derived by our numerical procedure. 

In the following, for simplicity we are assuming a thermal average for the cross-section times Moeller 
velocity, given by (crv) = a + bx, in which S-wave and/or P-wave annihilation processes are important. 
In all cases discussed the value of the constant c/ in Eq. (fl"8| is taken c/ = V2— 1. When dilaton effects 
are absent this value gives the best fit, better than 5%, matching very well numerical and approximate 
results. In the presence of the dilaton, better convergence of (fll)l) to the numerical result is obtained 
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if we use the freeze-out point x c j defined by the criterion n (a v) = H. This is actually implicit in the 
calculation that led to Eq. (TTB]). For values of Cf that are of order unity we have checked numerically 
that is smaller than Xf, derived from Eq. (|18[) . having a value x^ = Xf / ' s where the shift coefficient s is 
close to s = 1.3. Using the shifted value Xf, instead of Xf, in Eqs. (fl7|) yields results for the relic density 
that approximate better the numerical values for the relic density as we shall see. Note that the use of x c j 
in (|17[) . it enhances R by a factor of about two not spoiling therefore the exponential suppression caused 
by this factor to the relic density. 

Sample results are presented in Table U where in the first two columns we give the values of the 
coefficients a, 6, in cm 3 /s, pertinent to the thermal averaged cross-section, and in the third column the 
WIMP mass. In the next four columns we display in the following order, the "shifted" freezing point, 
x^, the reduction factor 10 2 x R, the relic density as this is derived by solving numerically Boltzmann 
equation, Qhy \ num , and its approximate value O/ig \ app . In the last column the ratio flhy \ num / flh^ \ app 
is shown. We warn the reader that the values of the relic densities appearing in the Table are only 
for demonstration and are not within the WMAP range. In deriving ilh^ \ app , the factor k appearing 
in Eq. ( ITT]) has been approximated by taking the average value £(x) — (£(#/) + £(xx))/2 within the 
integral J* 1 £(x) (av)dx which defines the denominator of k, as has been already discussed. Then 
k = 2 t;(x f) / (£,(x f) +£(xx)) which as we have verified, following our numerical procedure, approximates 
fairly well the numerical result for the same quantity. Notice that in this approximate form the factor 
k is a function of the ratio £{xx)/£{xf) which can be easily expressed, using Eq. ((6]), in terms of the 
parameters b,c, Ax and the freeze-out temperature Tf . The results displayed in Table U show that the 
semi-analytic scheme is successful in reproducing satisfactorily the numerical results with an accuracy in 
the range ~ 15%. Therefore although not as accurate as in the conventional case, where semi-analytic 
schemes yield better accuracies, the approximation used in the case at hand gives a fairly good picture 
of the relic density when dilaton effects are taken into account. 

Concerning the factor R a few remarks are in order. The values of R obtained are of the order of 
~ 0(1O -2 ) and for fixed input values of a,b decrease with increasing the mass m x . For high values of 
rn^, above a few TeV, R can become as small as ~ 0(1O~ 3 — 1CP 4 ). For small masses, smaller than 
50 GeV or so, the exponent in Eq. (jTTJ) may not be small enough, being overwhelmed by the presence 
of £~ l (xf) in Eq. (TT71) . This occurs because in these cases the point xx = Tx/m x , which sets the lower 
limit of the integration in the exponent of Eq. (fTT]), becomes larger getting closer to the upper limit of 
the same integration. In these cases therefore the factor R may be larger than unity causing no-reduction 
but enhancement instead to the relic density, by a factor which is ~ 0(2). 

For supersymmetric models, which we discuss in the following section, the reduction factor causes 
dilution by factor of order 0(5 — 50) , the smaller (larger) corresponding to lighter (heavier) neutralino 
masses. The amount of the dilution depends on the particular SUSY inputs and therefore it is important 
that we scan the supersymmetric parameter space to delineate regions compatible with the current 
accelerator and cosmological data, including those from direct dark matter searches. This we do in the 
following section. 

IV. PREDICTIONS OF THE CMSSM 

As it was explained in the previous section the dilation terms may have significant impact on the 
predicted value for the neutralino relic density. In general, as shown in the previous section, for a specific 
point of the parameter space, the neutralino DM abundance is reduced. The effect of the reduction factor 
is plotted in Fig. [51 where in the context of the CMSSM we present examples of the (mo, M-y n) parameter 
space for tan/3 = 5 and 10 respectively, taking the trilinear coupling A = 0. The purple thin lines, that 
they cross almost vertically the plane, mark the magnitude of this reduction factor in the relic density. 
Along the dotted purple line on the far left this factor is 0.1, along the solid line it is 0.06 and finally on 
the dashed line on the far right this is 0.006. 

Turning to other details in these figures, the thin yellow shaded region on left is excluded by the LEP2 
searches (chargino mass bound). The dark red region, of triangular shape, at the bottom is excluded 
because there the stau is the LSP. In the hatched grey regions at the upper left corner, whenever they 
appear, the EW symmetry breaking is not successful. The very thin blue strip that lies above the stau 
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FIG. 2: The (mo, M1/2) parameter space for tan/3 = 5 and 10, assuming Ao = 0. The lines and regions are as 
described in the main text. 



LSP region having also a segment that follows closely the edge of the EWSB forbidden region, delineates 
regions that are compatible with the WMAP7 Q data if dilaton effects are ignored. On the other hand, 
the dark green region is the WMAP7 compatible region in the presence of the dilaton, while within 
the light green region the neutralino relic density is smaller than this bound. One immediately notices 
the significance of the new dilaton terms in the Boltzmann equation, by the resulting enlargement of 
the cosmologically favored regions. In the figures displayed the blue dashed-dotted curve denotes the 
light Higgs bound from LEP2, mu = 114 GeV. The part of the parameter space on the left of the solid 
black line, marked by LHC, is already excluded by LHC experiments [6j. However large portions of the 
cosmologically allowed regions have been now moved to regions (dark green) of the parameter space that 
are not excluded by the current LHC data. 

Recent CERN-LHC data from Higgs searches provide additional constraints. The absence of pseu- 
doscalar decay A — > t + t~ in ATLAS and CMS experiments puts bounds on the pseudoscalar mass rriA 
for given tan/3 (46|. In particular, for large tan/3 ~ 50 values ttla < 450 GeV are excluded. However 
this bound is not relevant in constrained supersymmetric scenarios, like CMSSM and mSUGRA, due to 
the fact that the masses of the Higgs bosons, but the light one, are quite heavy unless the SUSY breaking 
parameters are confined in the region of very small 7%, M 1 / 2 values. Therefore, this bound is practically 
inoperative in our analysis. Recall that large tan/3 ~ 50 values allow rapid neutralino dark matter anni- 
hilations via the pseudoscalar A-resonance in the conventional mechanism where dilaton couplings to DM 
are neglected. This region has been phenomenologically studied recently, in view of the on-going LHC 
experiments, since the channel pp — > b + A (or H) — > b fi + fi~ + X allows for amplification of the signal to 
background-ratio providing unique possibilities for measuring the Higgs boson masses and decay widths 
47] at CERN LHC experiments. 

As for the light Higgs-boson, both LHC experiments have reported possible evidence of a Standard 
Model like Higgs in the mass range 123 - 128 GeV 0- Thus in the figures presented in this paper we 
plot the curve along which rrih — 123 GeV, bearing in mind that in the regions of the CMSSM that we 
are studying the light CP-even Higgs boson accurately mimics the Standard Model Higgs. In Fig. [2] for 
tan/3 = 5 and 10 respectively, all the displayed parameter space yields Higgs masses smaller than 123 
GeV. For larger values of tan/3, displayed in the other figures, we notice that large part of the parameter 
space is compatible with nth > 123 GeV, which lies in the region right to the blue curve that illustrates 
the mh = 123 GeV Higgs boson mass. Actually, in Fig. |3]presented here, with M 1 / 2 < 3 TeV, the Higgs 
boson is always lighter than 124-125 GeV. Note that a light Higgs boson with mass larger than 120 or 
123 GeV it is not easy to be compatible with cosmologically favored regions of CMSSM for all the values 
of tan /3, as was discussed recently j48|. But, taking into account the dilaton effects on the computation of 
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FIG. 3: Like the Fig. [2] for tan [3 = 30 and 40. 
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FIG. 4: Like the Fig. [2] for tan p 



and 50. 



the neutralino relic density, we notice that large parts of the parameter space become compatible both to 
WMAP bound and to Higgs mass larger that 123 GeV, especially for tan/3 > 25. In particular, for tan (3 
in the range 25 to 40 large portions of the "dilaton deformed" coannihilation corridor have Higgs masses 
larger than 123 GeV, as can be seen in Fig. [3] For tan /3 in the range 40 to 48 there are such regions that 
are the "dilaton deformation" of the pseudo scalar Higgs boson funnel, which is the blue shaded region as 
shown in Fig. [4] left. It is obvious, that larger regions are compatible to these Higgs boson masses if one 
assumes that neutralino relic density is smaller than the WMAP bound, which allows for contribution to 
dark matter density from other particle species. 

The cyan dashed line, on the left of each figure, delineates the bound set by XENON100 direct detection 
measurement [30| . The excluded region is on the left of this curve and includes the focus point region 
in cases where the dilaton effects are ignored. It's worth noticing that in order to apply the neutralino 
DM direct detection limit, one must take into account the hadronic uncertainties that dominate the 
neutralino-nucleon interaction. The dominant uncertainty stems from the strange quark density in the 
nucleon (AT|ss|iV), which reflects the uncertainty in the pion nucleon sigma term E^jv = l/2(m u + 
rrid) (N\ uu + dd \N) [49j . Various estimates of the S^at vary from 64 ± 8 MeV to the more conservative 



14 



50 ± 14 MeV [sol ] ■ In this analysis we have adopted the value S^at = 50 MeV. Even with this conservative 
value a significant portion of the focus point corridor, especially for large tan /?, is excluded. The reason is 
that large part of the WMAP7 cosmologically allowed region, in the absence of the dilaton, is within the 
XENON100 exclusion area, as is seen, for instance in Fig. [4] This happens for values of M 1 / 2 up to 1000 
GeV, that correspond to m x < 450 GeV, where the XENON bounds are fully applied. 2 . Interestingly 
enough, the dilaton effects reshape the WMAP allowed region in such a way that it now lies in a region 
where XENON100 bounds are evaded. Hence the focus point region is a viable possibility in this scenario. 

Concluding this section, we showed that the dilaton coupling to DM may dramatically modify the 
predictions of the supersymmetric models. In the regions of the CMSSM parameter space that is accessible 
to the LHC experiments up to the y/s = 14 TeV, the combination of WMAP7, the Higgs boson searches 
at LHC and the direct detection experiments like XENON100, put new phenomenological restrictions 
[48| . However in the presence of dilaton driven terms in Boltzmann equation these regions are deformed in 
such a way that much larger parts of the supersymmetric parameter space arc compatible to cosmological 
bounds, to XENON100 and to Higgs boson masses in the region 123-128 GeV, as it is indicated from 
the recent LHC data. This opens new directions for further phenomenological analyses. 



V. INDIRECT DARK MATTER SEARCHES: 7- RAY FLUXES 

Indirect DM searches, dedicated to detect potential 7, neutrino or positron fluxes, produced by the 
DM particle annihilations, have been intensified during last years. Especially, photon data collected 
from various regions of the Universe, like the centers of galactic halos where a large density of DM is 
clumped, if they are above the known backgrounds, may signal the presence of DM. Currently these 
can be detected by the new generation of detectors whose experimental precision has been increased 
significantly, like Fermi-LAT or HESS, and ground-based experiments expected to operate in the future, 
like the Cherenkov Telescope Array (CTA). 



and HESS [31j experiments, and 
which is scheduled to operate in 



In this part using the 7-ray data, delivered by Fermi-LAT [21], |51 
taking into account the projected sensitivity of the CTA detector 
the future, we briefly discuss the consequences of the dilaton scenario for the aforementioned indirect DM 
searches, and investigate how dilaton modifies the conventional picture. In the analysis of the previous 
section we have mainly concentrated on predictions of supersymmetric models, namely the CMSSM, and 
saw that the dilaton may dilute the relic density altering the constraints imposed by the WMAP data, 
allowing smaller neutralino annihilation rates. This in turn, implies that supersymmetric models are not 
expected to be tightly constrained by indirect Dark Matter search data. However given the fact that in 
the near future new more precise experiments will start delivering data, improving the limits put on the 
DM annihilation cross section, it is important that we consider current and proposed limits put by these 
experiments too. It should be stressed that such an analysis is important not only for supersymmetry, 
but for any model that predicts the existence of Cold Dark Matter that annihilates to ordinary matter, 
irrespectively of the nature of the DM particle and the particular mechanism through which annihilation 
into standard model particles proceeds. 

Upper bounds on the WIMP pair annihilation cross section can be obtained from the 7-ray flux, 
due to the WIMP annihilations in various regions of the galactic DM halo [3(1 [55Tj57j . This flux is 
proportional to the non-relativistic limit of the thermally averaged annihilation cross section times the 
relative velocity (eri>), the energy spectra dN/dE 1 originating from DM annihilation into final states 
qq,T + r~ , /i + /i~ ,W + W~ , ZZ, Zh etc, and the J-factor that depends on the DM profile. In particu- 
lar, the gamma-ray signals from dwarf spheroidal galaxies observed by Fermi-LAT [53[ exclude WIMPs 
with masses m x < 40 GeV annihilating into bb,T + r~ , while HESS observations provide constraints for 
heavier WIMPs [5a, [58[. In order to study the indirect DM searches, in conjunction with the dilaton 
scenario, we will concentrate on various 7-ray data originating either from the galactic center (GC) or 
the dwarfs spheroidal (dSph) galaxies. In Fig. [5] we present constraints from Fermi-LAT, HESS and 



2 The published XENON100 bound applies for m x up to 1 TeV, although an extrapolation of the bound can be performed. 
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FIG. 5: Bounds of various indirect gamma ray experiments on the ((a v) , m x ) plane. The profires EIN = Einasto 
and NFW = Navarro-Frenk- White are used. 

CTA (projected) on the the ((cv) ,m x ) plane, assuming that the dominant DM pair annihilation goes 
through the bb channel and all other channels have negligible contribution. We also assume that since 
the WIMPS annihilating in the galactic halo today are very non-relativistic, the annihilations are es- 
sentially pure S-wave. Both of these assumptions are realized in the major portions of the parametric 
space of supersymmetric models but evidently can hold in other cases too. The Fermi-LAT experiment 
can detect 7-rays in the energy range 20 MeV to 300 GeV. For this analysis we use the Fermi-LAT data 
recorded between August 4th, 2008 and April 20th, 2012. We perform a data analysis using the package 
ScienceTools as it is suggested by the collaboration. As it was done in [53[ we are using data from 
ten dSph galaxies: Bootes I, Carina, Coma Berenices, Draco, Fornax, Sculptor, Segue 1, Sextans, Ursa 
Major II and Ursa Minor. Reproducing the statistical analysis in (53|, and using the J- factors for dSph 
galaxies used in this analysis, we calculate the combined 95% CL upper limits shown in Fig . [5] for the 
joint likelihood of these ten dSph. These are in agreement with Fig. 2 presented in 53]. Looking 
into the physics of this constraint, one can see that practically assuming that the dominant annihilating 
channel is the b b only models with low m x are constrained. In fact, in the presence of the dilaton, 
models with m x up to 25 GeV are excluded. This conclusion is drawn using the WMAP compatible 
bands that are plotted in Fig. [5] Within the light blue horizontal band of the {{uv) , m x ) plane the relic 
density is within the WMAP limits in the conventional case, i.e. ignoring dilaton effects, and the green 
strip bending downwards as mass increases designates the WMAP compatible region when the dilaton is 
present, as was described in section III. Notice that for m x < 30 GeV the dilaton predictions requires 
larger cross sections than in the conventional scenario and the green strip lies above the corresponding 
blue one. In fact, in this low m x regime the factor R causes enhancement, as has been explained in the 
previous sections. 

In the same figure we have also plotted the regions that are excluded by HESS measurements, based 
on 7-ray fluxes originated from the GC [58j . assuming an observation window of 10 degrees around this. 
The HESS instrument is sensitive to photons in the energy range 300 GeV to 30 TeV. In order, to 
delineate the exclusion region at 95% CL on the ((erf) , m x ) plane, we divide the HESS energy range 
into 30 bins evenly spaced on a logarithmic scale. For the calculation of the supersymmetric signal 
we include the monochromatic, the continuum and inverse Compton scattering components as they are 
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FIG. 6: Like Fig. [4] including contours of constant (cry), in cm 3 /s. Also shown are the regions expected to be 
excluded by CTA experiment from the projected data collected from measurements around the GC, as described 
in the text. 



described in detail in [fuj. We use two profiles for the halo density around the GC: the Navarro- Frenk- 
White (NFW) [U and Einasto [H, with the choice p = 0.3 GeV/cm 3 , for the local DM density, and 
Ro = 8.5 Kpc for the solar distance to the GC (more details on the halo profile parameters and the 
methodology that we follow can be found in [60]). The background contribution consists of the diffuse 
galactic emission, the isotropic extragalactic contribution and resolved point source contribution. We 
model these background contribution using the latest models that are available from the Fermi-LAT 
collaboration [63|. We construct a \ 2 defined by \ 2 — J2i Sf /{Si + Bi), where Si and Bi are the signal 
and the background counts per bin, respectively. Then we compute the corresponding 95% CL exclusion 
limit based on this x 2 which gives a curve quite similar to the bounds presented in Fig. 4 of [58| . In the 
same manner we can draw the region on the {(crv) , m x ) space that can be reached by the future CTA 
detector (54|. We use a similar \ 2 analysis, but calculated with an observation time and effective area 
for this detector as given in Table II of [36|], following the discussion presented in section V of this paper. 
The projected excluded region by CTA is again presented in Fig. [5] for the NFW and the Einasto profiles. 
This is similar to Fig. 2 of [36[ when one uses the same units and divides the cross section by to 2 . Based 
on these, we are expecting that CTA detector can probe annihilations with values of (a v) of the order 
of 10 _26 cm 3 /s that overlap, in the conventional model, with the cosmologically favoured WMAP region 
for m x masses in the range m x ~ 400 GeV — 5 TeV as can be seen from Fig. [5] On the other hand, 
in the context of the dilaton cosmologies, discussed in previous sections, much smaller values of (a v) are 
compatible to the WMAP data and therefore in the presence of the dilaton studied in this work, it seems 
that even with the CTA scheduled performance it will be hard to derive a useful constraint. 

Passing to the case of CMSSM, in Fig. [B]we present in the (mo, Mj/j) parameter space, for tan j3 = 48 
(left plane) and 50 (right plane) , contours of the total pair annihilation cross section (a v) for values 
of it 10 -28 , 10 -27 , 10~ 26 and 5 x 10~ 26 cm 3 /s respectively. We notice, as expected from the previous 
discussion on Fig. [5j that the parts of the parameter space that are cosmologically favoured in the 
conventional model are compatible with values (crv) ~ 2.5 x 10 _26 cm 3 /s. Therefore these regions can 
be probed by the CTA expected measurements. However, these measurements will be unable to probe 
models with relic density calculations affected by the dilaton effects. To illustrate this in greater detail 
we also mark on these figures the regions of the parameter space likely to be constrained by the CTA 
from the GC (10 degrees window) for the NFW (cyan inclined lines) and the Einasto (blue horizontal 
dots) profiles. These regions surround the cosmologically allowed regions of the conventional model, 
but stay apart from the corresponding WMAP allowed regions in the presence of the dilaton, being in 
qualitative agreement with the findings regarding Fig. [5] Note also that in this case for the calculation of 
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this projected bound we have used all the neutralino pair annihilation channels occurring in the CMSSM 
model and not just the bb as we did before when considering the case depicted in Fig.[SJ The annihilation 
channel bb still dominates in the major portion of the parameter space, except the focus-point region, 
where W + W~ dominates, and the neutralino-stau coannihilation region where t + t~ is the dominant 
channel. 

The conclusion reached in this part is that the cosmologically favored regions of the parameter space 
of the CMSSM can probably beyond reach of future planned indirect 7-ray DM detection experiments, 
due to the dilaton effects which dilute the DM relic density preferring smaller annihilation cross sections. 
Although we have demonstrated this for high values of tan/3 this holds true for smaller values of the 
parameter tan/3 as well. In particular in these cases, the regions compatible with the cosmological 
bounds for the relic density have cross sections (crv) ~ 10~ 28 cm 3 /s, or less, well below the sensitivities 
of current and planned future experiments. 

VI. CONCLUSIONS 

In this paper we have shown that the dilaton dynamics during early eras, long before Nucleosynthesis, 
in conjunction with its coupling to Dark Matter may have dramatic consequences for the predicted Dark 
Matter relic density in the popular supersymmetric models. Modeling the dilaton evolution to be that 
dictated by exponential type potential, occurring in quintessence scenarios, inflation models and string 
theory, the ordinary predicted DM density may be diluted by large factors ranging from 0(5) to O(50) . 
This dilution mechanism is consistent with the absence of dilaton couplings to ordinary matter ( hadrons ), 
in the continuity equations, but it affects DM relics since dilaton dominates over radiation during and after 
DM decoupling. This allows for LSP annihilation cross sections, in the popular supersymmetric schemes, 
that are smaller by an order of magnitude or more. This however may imply smaller inelastic cross 
sections of the neutralino LSP with nucleons puttin g fa rther the potential of discovering supersymmetric 
DM at proposed direct detection experiments [29l. l30j] . We studied the predictions of the CMSSM in 
the light of LHC recent data and cosmological constraints arising from WMAP7 and XENON100 which 
puts the most stringent constraints of all direct DM search experiments. We found that the allowed 
points cover broad regions that occupy different regions of the parametric space as compared to ordinary 
models. This opens new directions for further phenomenological studies and prolongs the viability of 
supersymmetric models. 

As for indirect detection experiments, indirect searches of Dark Matter through antimatter produc- 
tion has stirred much interest the last three years. These data may be conditionally explained as DM 
annihilation in the galactic halo, generating the produced antiparticle flux. However more conventional 
interpretations exist that can explain the origin of these fluxes. The 7-rays measured by various experi- 
ments is another powerful tool for indirect detection of DM. In the cases studied in this work the smaller 
annihilation cross sections, which are required in this dilaton scenario in order to satisfy the WMAP 
data, weakens the potential of discovering neutralino DM in the galactig halo, through 7-ray current and 
future planned experiments, unless the detector performances are greatly improved. 
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